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Abstract 

We report the observation of ferromagnetic resonance (FMR) and the determination of the 
magnetocrystalline anisotropy in (lOO)-oriented single-crystalline thin film samples of Gai-^Mn^P 
with x = 0.042. The contributions to the magnetic anisotropy were determined by measuring the 
angular- and the temperature-dependencies of the FMR resonance fields and by superconduct- 
ing quantum interference device magnetometry. The largest contribution to the anisotropy is a 
uniaxial component perpendicular to the film plane; however, a negative contribution from cubic 
anisotropy is also found. Additional in-plane uniaxial components are observed at low tempera- 
tures, which lift the degeneracy between the in-plane [Oil] and [Oil] directions as well as between 
the in-plane [010] and [001] directions. Near T = 5 K, the easy magnetization axis is close to 
the in-plane [011] direction. All anisotropy parameters decrease with increasing temperature and 
disappear above the Curie temperature Tq- A consistent picture of the magnetic anisotropy of 
ferromagnetic Gai-^Mn^P emerges from the FMR and magnetometry data. The latter can be 
successfully modeled when both coherent magnetization rotation and magnetic domain nucleation 
are considered. 

PACS numbers: 75.30.Gw, 75.50.Pp, 76.50.+g 
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I. INTRODUCTION 



Mn-based diluted magnetic semiconductors show a variety of different magnetic ordering 
phenomena ranging from ferromagnetism mediated by quasi-delocalized holes in materials 
exhibiting metallic conductivity [l|, |2[|_to spin-glass like behavior in semiconducting matrices 
attributed to Mn-rich nanoclusters ^j]. An important parameter expected to govern the 
magnetic ordering is the localization of the charge carriers coupling the 3d high-spin states 
commonly introduced by Mn incorporation. In III-V materials, where Mn simultaneously 
acts as acceptor, the corresponding acceptor level essentially determines the degree of local- 
ization of the holes \&\. A variation of the acceptor level can be achieved by changing the 
group-V atom in III-V alloys. While recent studies on Gai^Mn^N indicate the formation 
of a ferromagnetic ordering in essentially insulating material with a Curie temperature Tq 



of 8 K 



5j , a carrier- mediated, non- metallic phase with Tq up to 65 K in Gai_a;Mn x P has 



recently been synthesized [6 




i El- 



the magnetic dopant concentration and that ferromagnetism is suppressed by the addition 
of compensating Te 6] and S donors jg|. Furthermore, x-ray absorption spectroscopy and 
x-ray magnetic circular dichroism have indicated that the local ferromagnetic environment 
for Mn in Gai-zMn^P is very similar to that in Gai-^Mn^As and that the Mn <i-derived 
density of states at E F is strongly spin polarized {{J. 

In Gai_ x Mn x As, the magnetocrystalline anisotropy has been successfully described in 
terms of the GaAs valence band because the states occupied by holes mediating inter-Mn 



exchange appear to be sufficiently simi 



has been previously pointed out 
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ar to the unperturbed GaAs valence band 10[. It 
12j that the holes responsible for exchange mediation 
are probably at least semi-localized in real space, as assumed in polaronic theories Q] 
also used to describe diluted magnetic semiconductors. Evidence for some degree of hole 
localization in III-Mn-V ferromagnetic semiconductors has been observed in infrared stud- 
ies of low-temperature molecular beam epitaxy (LT-MBE) grown Ini_ x Mn x As 15j and 
Gai-^Mn^As 16J. In these materials it is believed that there is significant mixing be- 
tween the valence and Mn impurity bands, while our studies in Gai^^n^-P suggest that 
the exchange- mediating holes reside in a separate Mn impurity band jy, 7, sj]. Conduction in 
this band occurs by hopping in the ferromagnetic regime, indicating a much higher degree of 
localization than in Gai_ a Mn x As or Ini.^Mn^As [zj]. This suggests that, as predicted by ab- 
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initio calculations 
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211 ] , a discontinuous transition to an alternative mechanism 



of carrier-mediated exchange does not occur even with these distinctions in band structure. 
As the magnetic anisotropy is intimately tied to the properties of the exchange-mediating 
holes, it is important to investigate the magnetic anisotropy connected with the impurity 
band in Gai-^Mn^P. Therefore, after giving a short introduction into the sample fabrica- 
tion and the measurement techniques in Sec. [Til we determine the contributions to the mag- 
netic anisotropy of Gai_ x Mn x P by measuring the angular- and the temperature-dependence 
of ferromagnetic resonance (FMR) fields in Sec. IIII1 The FMR results are substantiated 
by superconducting quantum interference device (SQUID) magnetization measurements in 
Sec. EH 



II. EXPERIMENTAL 



Samples were prepared by ion implantation followed by pulsed-laser melting (II-PLM) 22, 



231 ] . For the present study, unintentionally sulfur-doped n-type GaP (100) wafers with a 
carrier concentration in the range of 10 16 cm" 3 to 10 17 cm -3 were implanted with 50 keV Mn 
ions. Each implanted sample was cleaved to have [Oil] and [Oil] edges and was irradiated in 
air with a single 0.4 Jem -2 pulse (FWHM = 23 ns) from a KrF excimer laser (A = 248 nm) 
homogenized to a spatial uniformity of ±5% by a crossed-cylindrical lens homogenizer. 
Channeling 4 He + Rutherford backscattering spectrometry (RBS) and particle induced x-ray 
emission (PIXE) were used to assess the crystalline quality, Mn dose retained after II-PLM, 



and substitutional fraction of Mn in the samples [23j. Once processed the films are high- 
quality single crystals with a Mn dose of 7.3 x 10 15 cm -2 and substitutionality of 0.7, i.e., 70% 
of Mn atoms substitute Ga atoms. We note that the remaing 30% of Mn atoms do not form 
interstitial defects and instead are incommensurate with the GaAs lattice presumably in 
the form of small clusters. This level of substitutionality is not unlike Gai-xMn^As films of 
higher Mn concentration grown by low-temperature molecular beam epitaxy, which contain 
on the order of 20% non-substitutional Mn 241 ] . 



II-PLM processing results in samples having a gradient in Mn concentration into the 
depth of the sample as measured by secondary ion mass spectrometry (SIMS), making it 
impossible to determine single values for the film thickness and Mn concentration. The Mn 
SIMS profile can be approximated by a Gaussian distribution centered at a depth of 40 nm 
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FIG. 1: Temperature dependence of the FMR signal of a Gao.95sMno.042P sample for the magnetic 
field aligned along the in-plane [Oil] (left panel) and the out-of-plane [100] (right panel) directions. 
The dashed lines are guides to the eye. The magnetic field corresponding to g = 2 is indicated by 
dotted vertical lines. 

with a width of 20 nm. However, as the regions of the film with highest Mn concentra- 
tion dominate both the magnetic and transport properties, samples are discussed here in 
terms of their peak substitutional Mn concentration, which was determined by channeling 
RBS and PIXE to be x = 0.042 Sample magnetization was determined in various 

crystallographic orientations using a superconducting quantum interference device (SQUID) 
magnetometer. The FMR measurements were performed at iojlix « 9.3 GHz in an elec- 
tron paramagnetic resonance (EPR) spectrometer using magnetic field modulation, with the 
sample temperature controlled using a liquid-He flow cryostat. 

III. FERROMAGNETIC RESONANCE SPECTROSCOPY 

The left panel of Fig. [1] shows the temperature dependence of the FMR signal of a typical 
x = 0.042 sample for the magnetic field aligned along the in-plane [011] direction, while 
the right panel presents data from the out-of-plane [100] direction. At T = 5 K for H 
parallel to the in-plane [011] direction, we observe one resonance at fioH = 236 mT with 
a peak-to-peak linewidth of fioAH pp ps 36 mT, while for H perpendicular to the sample 
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FIG. 2: Angular dependence of the ferromagnetic resonance fields for the magnetic field rotating 
within the (Oil) (left panels) and (100) (right panels) planes at T = 5 K, 15 K, 25 K, 35 K, and 
45 K. The circles correspond to the experimentally observed resonance positions. The full lines 
show the anisotropy expected for the parameters shown in Fig. UJ 
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FIG. 3: Coordinate system used for FMR simulation, with the orientation of the saturation mag- 
netization M = M(6, <£) and the magnetic field H = H(6, <p). 

surface (H || [100]) there are three distinct resonances at fioH = 576 mT, 492 mT, and 
430 mT, each with /!oAif pp m 25 mT. With increasing temperature, the resonance fields 
for both orientations shift toward fioH = 330 mT, which corresponds to the resonance field 
of paramagnetic impurities with a ^-factor of g = 2. The anisotropy disappears around 
T pa 65 K, slightly above the Curie temperature Tq = 55 K determined from SQUID 
magnetization measurements. We attribute this to the moderately-large applied field hqH 
330 mT in resonance, which stabilizes ferromagnetism even slightly above Tq (compare Fig. 3 
in Ref. [9|). 

We attribute the multiple resonances in the H || [100] data to spin wave excitations 



25 



26 



27 



. Since we only observe three resonances or less, a detailed analysis of the mode 
spacing is hardly possible. However, assuming the resonance at the highest magnetic field 
IIqHq = 576 mT to be the collective mode, we determine a separation S(n) = Hq — H n oc 



n 



0.8 



between the resonance fields of the mode with the highest field Hq and the nth spin 



wave mode H n at T = 5 K, which does not obey the classical behavior expected for a 
homogeneous film E c i asgica i(n) oc n 2 [26|. We attribute this to the varying depth profile of the 



Mn concentration. A similar non-quadratic behavior of the mode spacin g ha s been reported 



25 



27 



28 



. Describing 



for Gai_ x Mn x As thin films exhibiting gradients in hole concentration 
the implantation profile with a parabolic depth dependence of the Mn concentration, a linear 
modes spacing would be expected, in reasonable agreement with the observed behavior. 

To elucidate the magnetic anisotropy of the samples we performed measurements of the 
angular dependence of the FMR for sample rotations in the (011) and (100) planes. The 
angular dependences of the resonance fields obtained at different temperatures are shown in 
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Fig. [2j The panels on the left hand side correspond to rotations of the external magnetic field 
from the in-plane [Oil] to the out-of-plane [100], and back to the in-plane [Oil] direction. 
For simplicity, we limit the discussion to the collective mode in the following. A uniaxial 
magnetic anisotropy with the magnetic hard axis perpendicular to the layer can be inferred 
from the increase in ferromagnetic resonance fields approaching [100]. Likewise, the fourfold 
symmetry observed for the in-plane rotations in the right hand panels demonstrates the 
presence of a cubic anisotropy contribution. Additionally, the fact that the [Oil] and [Oil] 
orientations, as well as the [010] and [001] orientations are not degenerate indicates the 
contribution of further in-plane uniaxial anisotropy components. 

For a quantitative simulation of these data we use the free energy density 

F— — Mif (sin sin $ sin 9 sin <fi + cos G cos 6* + sin G cos $ sin# cos0) 
+ Klg° sin 2 G sin 2 $ 

- -K^x sin 4 G sin 4 $ 

- -K^ (cos 4 G + sin 4 G cos 4 $) 
+ ii^ n (cosQ + sinQcos<!>) 2 

+ K° 01 sin 2 Gcos 2 $. (1) 

The angles are given by the orientation of the saturation magnetization M = M(Q, $) 
and the magnetic field H = H(6,<f)) (Fig. [3]). The first term describes the Zeeman energy, 
while the second term represents an effective perpendicular uniaxial anisotropy and 
is composed of the sum of demagnetization and uniaxial magneto crystalline components, 
|/ioM 2 and K^ 00 , respectively. In order to describe the breaking of the cubic symmetry due 
to in-plane biaxial compressive strain induced by the presence of Mn, we include separate 
cubic terms for in-plane and perpendicular components, given by and K^_, respectively. 
The in-plane cubic symmetry breaking is accounted for by the final two terms representing 
uniaxial contributions along [011] {K® 11 ) and [001] (K® 01 ). In all cases the appropriate 
anisotropy field is given by the ratio 2K/M. The come about of the different anisotropy 
terms, as well as of the equivalence of first order cubic and second order uniaxial anisotropy, 
are discussed in Appendix A. 



S 



Following the approach of Smit et a/. .|29l. l30j we obtain the equation of motion 

2 ~> r/q2\/Q2\ / a a \ 2~ 



M 2 sin 2 6 



<9$ 2 \de 2 J \d$de 



(2) 

*0,©0 



with the gyromagnetic ratio 7 = which has to be evaluated at the equilibrium orienta- 
tion of the saturation magnetization determined from 

^ F l-f=*o = ^ F le=eo = 0. (3) 

The solution of these equations yields the FMR resonance condition. The full lines in 
Fig. [2] are simulations of the measured data with anisotropy fields plotted in Fig. HI At 
5 K, the results for four samples with x = 0.042 are well reproducible, with 0.16 T < 
2K^°/M < 0.19 T, -0.10 T < 2K^/M < -0.06 T, -0.04 T < 2K\/M < -0.03 T, 
T < 2K° U /M < 0.012 T, and T < 2A'° 01 /Af < 0.02 T. The magnetic anisotropy of 
these films is clearly dominated by the uniaxial and cubic contributions perpendicular to the 
film. The in-plane cubic and uniaxial anisotropics along [011] and [001] are all somewhat 
smaller. Figure |4] depicts the decrease of all of the anisotropy components with increasing 
temperature and demonstrates that all components disappear above the Tq of 55 K of the 
films as expected. As discussed above, only a small 2K^° is observed even above Tc at 
65 K. 

From SQUID magnetization measurements we estimate the saturation magnetization in 
the most heavily doped part of the film — i. e. near the peak of the Mn distribution 
— to M = 37 kA/m. From this, we obtain an upper limit for the demagnetization field 
HqM of 0.05 T. Therefore, the demagnetization field constitutes only about one fourth of 
the effective uniaxial magnetic anisotropy field along [100], 2K^°/M = 0.19 T, determined 
from the simulation of the angular dependence of FMR. This strongly indicates the presence 
of a tetragonal distortion of the Gao.95sMno.042P layer after pulsed-laser melting causing 
the dominating contribution A* 00 as in the case of Gai_ x Mn x As, where a strong uniaxial 
magnetic anisotropy in the growth direction is commonly observed and is attributed to the 



31 



32 



33]. 



tetragonal distortion of the Gai-ajMn^As layer due to lattice-matched growth 
However, a quantitative detection of the distortion in case of Gao.95sMno.042P via x-ray 
diffraction turns out to be difficult due to the inhomogeneous Mn profile after pulsed-laser 
melting and the accompanied broadening of the diffraction peaks. 
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FIG. 4: Temperature dependence of the anisotropy parameters and the ^-factor obtained from the 
simulation of the angular dependence of the ferromagnetic resonance (full lines in Fig. [2|). 

One very interesting finding is that the magnetic easy axes of the in-plane cubic magnetic 
anisotropy are along [011] and [011] directions as opposed to the [010] and [001] directions 
commonly observed for LT-MBE grown Gax_3.Mn3.As. This observation gives rise to the 
negative sign of K\ in the Gao.95sMno.042P samples studied (K^_ is also negative). To 
the best of our knowledge, a negative cubic anisotropy has so far only been reported for 
Ini_3;Mna;As 34], [35|. Within the different models for carrier- mediated ferromagnetism in 



diluted magnetic semiconductors 1(3, |36|], the sign of the cubic anisotropy is predicted to 



oscillate with varying hole concentration. However, the applicability of these models for 
the strongly localized, impurity-band-like character expected for the holes in Gai^Mna-P 
remains an open question. Taking into account all anisotropy contributions, the global 
magnetic easy axis at 5 K is oriented close to the [Oil] direction. 

Interestingly, the ^-factor does not deviate significantly from g = 2. In Gai_a;Mn~As, 
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g was found to be an effective g-factor taking into account both the contributions of the 
Mn atoms and the hole subsystem [3?]]. Depending on the hole concentration p, Liu et al. 
observed a q-factor at 4.2 K between q = 1.80 for a sample with p = 1.64 x 10 20 cm -3 and 

n 

g = 1.95 for a sample with p = 1.24 x 10 20 cm 3 [37j]. The fact that g is found to vary 
from 1.95 to 2 for increasing temperature from 5 K to 65 K indicates that there is only 
an even smaller contribution of the hole subsystem to the effective (/-factor in the case of 
Gai_ x .Mn x .P. This is consistent with the observations and calculations in Ref. [?J indicating 
a small hole concentration of up to 10 20 cm -3 . 



IV. SQUID MAGNETIZATION MEASUREMENTS 

To substantiate the results obtained from FMR in the preceding section we also per- 
formed SQUID magnetization measurements. Figure [5] compares the M(H) magnetization 
curves obtained at T = 5 K for the external magnetic field oriented along the out-of-plane 
[100] (solid triangles) and the in-plane [Oil] (solid squares) and [Oil] (open circles) crys- 
tallographic axes. The square-like M(H) curve obtained for H\\ [Oil] also indicates that 
[Oil] is the easy magnetic axis at 5 K. Similarly, the large field of ~ 0.2 T required to align 
M along the magnetically hard [100] direction is due to the large out-of-plane uniaxial and 
cubic anisotropy contributions. In the following, we use the free energy ansatz of (CD) to sim- 
ulate these M(H) curves and especially to explain the kink observed for the magnetization 
measurement along [011]. 

To begin with, we focus on the M(H) curve for H\ | [100] (out of plane) shown in Fig. E(a) 
on a larger field scale. The dotted line is the curve simulated as discussed below for which 
we obtained the best agreement with the SQUID measurement using the anisotropy fields 
2K^°/M = 0.1 T, 2K&/M = -0.12 T, 2K\/M = -0.04 T, 2K° n /M = 0.005 T, and 
2K® 01 /M = 0.004 T in (CQ). For H\ \ [100] the simulated curve is predominantly determined 
by 2K^°/M and 2K^/M. Both these parameters agree with the ones determined from 
FMR, to within a factor of two which can be understood as follows. The FMR spectra for 
H oriented along the hard magnetic out-of-plane [100] axis features spin wave excitations. 
This not only indicates that there is a gradient in magnetic properties as already discussed 
above, but also that these modes are located at the region of the highest uniaxial anisotropy 
field 2K^°/M and therefore only probe the magnetic properties of this specific region. 
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FIG. 5: M(H) SQUID magnetization curves along the out-of-plane [100] (solid triangles) and the in- 
plane [011] (solid squares) and [011] (open circles) crystallographic axes. The lines are guides to the 
eye. The right and left vertical axes give the magnetic moment per substitutional Mn atom myin and 
the magnetization M of the sample in the region of highest Mn concentration, respectively. These 
values were deduced from the measured total magnetic moment m tot as described in Appendix B. 

SQUID magnetization measurements in contrast integrate over the magnetic properties of 
the whole layer. In this respect the agreement of FMR and SQUID within a factor of two is 
quite reasonable. 

The curvature of the simulated magnetization curve can be understood looking at the 
dependence of the free energy per Mn atom on the orientation of the magnetization in the 
(011) plane for different magnetic fields applied along [100] [Fig. E^b)]. The position of the 
solid circle corresponds to the equilibrium orientation of magnetization in the minimum of 
the free energy surface. For high magnetic fields H\ \ [100], the magnetization is also in the 
[100] direction, since then the Zeeman term is the dominant contribution in ([I]). This is the 
case for noH = 0.3 T in Fig. [6](b). For decreasing magnetic fields the magnetocrystalline 
anisotropy becomes increasingly important. This leads to the migration of the minimum in 
the free energy surface — and therefore also equilibrium orientation of the magnetization 
— in the direction of the [011] axis, which is the magnetic easy axis for zero magnetic 
field. The application of a magnetic field in the opposite direction (i.e. H\ \ [100]) in turn 
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FIG. 6: (a) Comparison of the M{H) SQUID magnetization curve measured for H along the out- 
of-plane [100] direction (full triangles) with the curve simulated using the anisotropy parameters 
given in the figure (dotted line), (b) Free energy per Mn atom as a function of the orientation of 
magnetization in the (Oil) plane for different magnetic fields applied along [100] calculated from 
(PQ) . The position of the solid circle corresponds to the equilibrium orientation of magnetization in 
the minimum of the free energy surface. For clarity, the curves are shifted vertically. 

tilts the magnetization more and more in this direction [Fig. EJ^b)]. Having determined the 
orientation of the magnetization depending on the field strength H, we obtain the simulated 
curve in Fig. E^a) via calculating the projection of the magnetization along the direction of 
the external magnetic field, which is the quantity measured by the SQUID magnetometer. 
The process of magnetization reversal described in this paragraph is called coherent spin 



rotation |38| 



For the simulation of the in-plane M(H) curves in Fig. El^a), in addition to coherent 
spin rotation, the process of non-coherent spin switching has to be considered following the 
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FIG. 7: (a) Simulated magnetization curves using the anisotropy parameters given in the figure. 
The solid, dashed, and dotted curves correspond to H\ \ [011], H\ \ [011], and H\ \ [100], respectively, 
(b) Free energy per Mn atom as a function of the orientation of magnetization in the (100) film 
plane for different magnetic fields applied along [011] and (c) along [011]. The position of the 
solid circle corresponds to the equilibrium orientation of magnetization in the minimum of the free 
energy surface. For clarity the curves are shifted vertically. 
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discussion of Ref. [38|. For H\ \ [Oil], spin switching is visualized in Fig. [7(b), where the 
dependence of free energy per Mn atom on the orientation of the magnetization in the (100) 
film plane is plotted for different magnetic fields. Decreasing the external magnetic field 
H\ \ [Oil] from its maximum value of 7 T to zero, the magnetization remains "trapped" in 
the global minimum at [Oil]. For negative fields, [Oil] turns into a local minimum, while 
[Oil] becomes the global minimum in free energy. Since the thermal energy fc#T at T = 5 K 
of 0.43 meV is three orders of magnitude larger than the energy needed to overcome the 
barrier between the two minima, there will be always some magnetic moments oriented 
along the direction of the global minimum. However, for the generation of a new magnetic 
domain with a magnetization along [Oil], the domain walls of the nucleus of this domain 
first have to be formed. The formation of these domain walls is energetically unfavorable, 
since the magnetic moments in the walls are oriented along the magnetic harder axes of the 
energy barriers. Therefore, magnetization reversal takes place only if the energy gain from 
tilting the magnetization into the direction of the global minimum accounts for the energy 
needed for the formation the domain walls of this new magnetic domain. To obtain the 
experimentally observed switching field of -2.5 mT found for M(H), with H\ \ [Oil], we have 
to assume a domain wall formation energy of AE 011 = 8.6 x 10~ 4 meV per Mn atom. 

The situation for H\ \ [Oil] is plotted in Fig. UJc). Decreasing the external magnetic 
field if 1 1 [Oil] from its maximum value of 7 T to zero, the magnetization first also remains 
"trapped" in the global minimum at [Oil]. However, approaching /j,qH = mT, [Oil] 
becomes the global minimum due to the presence of the uniaxial anisotropy field along 
[Oil]. Therefore, there will be a first switching into the [Oil] direction at positive fields and 
a second switching into the [Oil] direction at large enough negative magnetic fields. In the 
model, we assumed the same energy barrier AE 011 for both switching processes. Then, the 
parameters predominantly determining the switching fields are AE 011 and 2K® ll /M. For 
increasing AE 011 both switching processes occur later (at lower fields), while for increasing 
2K° n /M the first switching process takes place earlier (at higher field) and the second one 
later (at lower field). The best agreement with the experimentally observed switching fields, 
fioH = 1.5 mT for the first and [1qH = —3.5 mT for the second switching [compare Fig. [8], 
we obtained for the domain wall formation energy of AE 011 = 1.6 x 10~ 4 meV per Mn 
atom and the uniaxial anisotropy field 2K® ll /M = 5 mT. Note that in contrast to the 180° 
domain walls that must be formed in case of the magnetization reversal for H\ \ [Oil], the 
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FIG. 8: Temperature dependence of magnetization for the sample cooled down to 5 K in a field 
of (ioH coo i = 7 T and measured during warm up in a field of fj,QH measure = 1 mT along the same 
crystallographic axis as in Fig. [SJ 

domain walls that must be formed here comprise a 90° rotation of magnetization. Therefore 
it is not surprising that AE 011 < AE 0U . Typical values we obtain for the height of the 
energy barrier separating two local energy minima AF = 7 x 10~ 4 meV per Mn ion are 



one order of magnitude smaller than the value obtained 
reversal process in case of Ga-|_, T Mn T As by Liu et al. 
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wall model as discussed by Liu et al. 



or a perpendicular magnetization- 
|. Applying Kittel's Bloch domain 



to our in-plane magnetization reversal process 
for H\ \ [011] we obtain a size of a domain nucleus of approximately 5 /im. However, note 
that the assumption of a Bloch domain wall may not be justified in this case of an in-plane 
magnetized ferromagnetic film. 

Consequently, the kink in the M(H) curve for H\ \ [011] can be explained by the fact that 
due to the presence of the uniaxial anisotropy field along [011], [011] is not the global easy 
magnetic axis. At low fields a non-coherent spin switching into the global easy magnetic 
axis along [011] takes place, which causes a vanishing projection of M on the [011] direction. 
The fact that there still is a finite projection in Fig. [7(a) is caused by the presence of the 
uniaxial in-plane anisotropy field 2K® 01 /M which slightly changes the 90° angle between the 
two minima of the cubic anisotropy near [011] and [011]. 

It should be noted that our simulation only accounts for hysteresis effects caused by 
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the non-coherent spin switching described above; the additional hysteresis effects observed 
for example for the M(H) curve measured for H along the out-of-plane [100] direction 
are presumably caused by the pinning and de-pinning of domain walls at crystal defects, 
which is not included in our model. Furthermore, the magnetization is by far not saturated 
at magnetic fields of HqH = 0.02 T. Therefore, the saturation magnetization used in our 
model M = 26 kA/m is lower than the real saturation magnetization of M = 37 kA/m 
at fioH = 7 T. Internal stresses, or shape irregularities could explain the rounding of the 
magnetization curve at high fields 39|. However, in spite of the simplicity of the model, 
the M(H) magnetization curves along several crystallographic directions can be explained 
at least semi-quantitatively by the presence of the anisotropy fields determined from FMR. 

Finally, we discuss the temperature dependence of magnetization along different crystal- 
lographic orientations (Fig. [8]). In these measurements the sample is cooled down in the 
maximum available field HqH coo \ = 7 T. At 5 K, the field is switched to ^oH measure = 1 mT 
and the projection of magnetization along the field direction is measured during warm up of 
the sample. For H\ \ [Oil] (closed squares), we obtained the highest value for this projection 
at all temperatures in agreement with [Oil] being the easy magnetic axis. Accordingly, the 
projection along [100] (closed triangles) is very small, since [100] is the hard magnetic axis. 
For if 1 1 [011] (open circles) and T < 50 K, the projection lies in between the values for the 
preceding orientations. This can be explained by the fact that [011] is not the global easy 
magnetic axis in this temperature range. Above 50 K, the curve for H\ \ [011] overlaps with 
the one for H\ \ [Oil], which is in good agreement with the disappearance of the uniaxial 
anisotropy field 2K® U /M in this temperature range (see Fig. HJ), which lifts the degeneracy 
between the easy axes [Oil] and [011]. Furthermore, the Curie temperature Tq = 55 K 
deduced from the M(T) curves in Fig. [8] again nicely agrees with the temperature around 
which the anisotropy fields vanish in Fig. HI 



V. CONCLUSIONS 



In conclusion, we have investigated the field and temperature dependencies of the mag- 
netic anisotropy of Gao.95sMno.042P thin films synthesized by ion-implantation and pulsed 
laser melting using measurements of the angular dependence of both ferromagnetic reso- 
nance and SQUID magnetometry. The results of FMR and SQUID measurements including 
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coherent spin rotation and non-coherent spin switching can be understood quantitatively us- 
ing a relatively simple free energy model. Similar to Gai_ x Mn x As thin films, the magnetic 
anisotropy is dominated by a strong out-of-plane uniaxial contribution. Since the demagne- 
tization field can only account for about one fourth of this out-of-plane uniaxial anisotropy 
field, its most probable origin is in-plane biaxial compressive strain, which is also the case 
for Gai-zMn^As thin films grown epitaxially on a GaAs substrate. We also observe a cubic 
anisotropy contribution. However, the sign of this cubic anisotropy term is opposite to the 
one commonly observed for Gai-^Mn^As. While the latter finding could still be in agreement 
with Dietl's theory of hole-mediated ferromagnetism considering the significantly reduced 
hole concentration in Gai-^Mn^P compared to Gai-^Mn^As, it remains to be demonstrated 
that this theory can indeed be applied to material systems exhibiting highly localized holes 
such as Gai_ x Mn x P. Nevertheless, it is an very interesting observation that in spite of the 
highly localized character of the holes in Gai-zMn^P, all the magnetic properties (satura- 
tion magnetization, absolute values of the anisotropy fields, Curie temperature) are similar 
to those typically observed in Gai-^Mn^As. This constitutes an important constraint for 
theories attempting to explain carrier-mediated ferromagnetism in highly localized material 
systems. Finally, the observation of symmetry-breaking in-plane uniaxial anisotropy com- 
ponents similar to that seen in Gai-^Mn^As — where the origin still is under debate — 
indicates an intrinsic origin related to the hole-mediated ferromagnetic phase in III-Mn-V 
ferromagnetic semiconductors. Moreover, because our samples were fabricated via II-PLM 
— a form of liquid-phase epitaxy — explanations brought forward that invoke vapor phase 
growth processes including effects of surface reconstruction can be excluded as the origin of 
this in-plane symmetry breaking. 

ACKNOWLEDGMENTS 

We thank E. E. Haller for use of the ion implantation facilities and I. D. Sharp and 
J. W. Beeman for experimental assistance. The work at the Walter Schottky Institut was 
supported by Deutsche Forschungsgemeinschaft through SFB 631 and the Bavaria California 
Technology Center; the work at Berkeley by the Director, Office of Science, Office of Basic 
Energy Sciences, Division of Materials Sciences and Engineering, of the U.S. Department of 
Energy under Contract No. DEAC02- 05CH11231 and previous contracts. 



18 



APPENDIX A: ANISOTROPY ENERGY 



Magnetic anisotropy can e.g. be caused by dipole-dipole interaction, crystal fields, and 
spin-orbit coupling. Furthermore, uniaxial and biaxial strain can be an origin for magnetic 
anisotropies. While to date no comprehensive ab-initio understanding of magnetic anisotropy 
has been established, magnetic anisotropy can be efficiently described mathematically with 
the help of symmetry considerations. According to Chikazumi 40(, the free energy of a 
uniaxial anisotropy can be expressed by expanding it in a series of powers of sin 2 



F u = K ul sm 2 $ + K u2 sm 4 $ + 



(4) 



with the first and second order constants K u ± and K u2 , respectively, and the angle i? between 
the orientation of magnetization m = || and the anisotropy axis u. This can be rewritten 
using sin 2 $ = 1 — cos 2 so that 



F n = K ul (l - cos 2 ■&) + K n2 (l - cos 2 d) 2 + . . . 

K ul + K u2 ) + (-Km - 2K u2 ) cos 2 + K u2 cos 4 + 



(5) 



const. + K u i cos 2 + K u2 cos 4 + 



M_ 

M 



with K u i := —K u i — 2K u2 and K u2 := K u2 . Using m 

cosines of the cartesian axes a x = sin sin $, a y = cos G, and a z = sin cos $ from Fig. [3l 









Oiy 




{a, ) 




= sin 


( 





with the direction 



the first order uniaxial anisotropy contribution along u = 
obtained via 



1 



in ([!]) for example is 



F u 011 = (muf = ^Kl l \cos 6 + sin O cos $) 2 . 



(6) 



The free energy for cubic magneto-crystalline anisotropy in cubic systems according to 
Chikazumi 40| is given by 

F c = K c i (a 2 x a 2 y + a 2 a 2 z + a 2 z al) + K c2 a 2 a 2 a 2 z + . . . (7) 

With the addition theorem for the direction cosines 

<4 + a t + <£ = 1 - 2 ( a l a l + a l a l + a l a l) > ( 8 ) 
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([7]) can be transferred to 



F c = const. — 2^ci («i + oiy + o^) + K c2 a 2 x a 2 a 2 + 



(9) 



Therefore, the af terms in ([9]) link the first-order cubic anisotropy and the second-order 
uniaxial anisotropy given in (jHJ). Consequently, the latter formulation for first order cubic 
anisotropy is equivalent to a linear combination of second order uniaxial anisotropies along 
the cartesian axes, F^ 2 = K^ 2 af, i G {x, y,z}. 

In ([I]) we accounted for the tetragonal crystal symmetry via distinguishing in-plane and 
out-of-plane cubic anisotropies. Following the above discussion one could equivalently use 
a combination of a first order cubic anisotropy and a second order uniaxial anisotropy 
perpendicular to the film plane 



with the first order cubic anisotropy constant K c i = and the second order uniaxial 
anisotropy constant perpendicular to the film plane = — ~~ -^li)- 
Note also that because of 



only two of the three first order uniaxial anisotropy constants K^°, K®\°, and K®® 1 are inde- 
pendent. A first order uniaxial anisotropy can always be expressed by two other first order 
uniaxial anisotropies perpendicular to each other. Analogously, the second-order uniaxial 
anisotropy constants also are not independent because of Eq. flSJ). 

APPENDIX B: DETERMINATION OF MAGNETIZATION 

Due to the Mn implantation profile, the magnetization can not be calculated as usual 
via dividing the total magnetic moment m tot measured e.g. via SQUID magnetometry by 
the sample volume. In a first step we calculate the magnetic moment per substitutional Mn 
atom via 




(10) 



a 2 x + a 2 + a 2 z = 1, 



(11) 



(12) 
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where -DMn,retained = 7.3 x 10 15 cm -2 is the Mn implantation dose retained after II-PLM, A is 
the sample area, and f su bst = 0.7 is the fraction of substitutional Mn atoms derived via RBS 
and PIXE. To obtain an estimate for the magnetization of the sample in the region of highest 
Mn concentration we multiply the average magnetic moment per Mn atom with the peak 
substitutional Mn concentration x and the density of Ga lattice sites [Ga] = 2.47 x 10 22 cm -3 

M = m Mn ■ x ■ [Ga] . (13) 
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